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Wavelets for Elliptical Waveguide Problems 

M. M. S. Lira* H. M. de Oliveira* R. J. Ciiitra* R. M. Campello de Souza§ 


Abstract 

New elliptic cylindrical wavelets are introduced, which exploit the relationship between analysing filters and Floquet’s 
solution of Matliieu differential equations. It is shown that the transfer function of both multiresolution filters is related 
to the solution of a Mathieu equation of odd characteristic exponent. The number of notches of these analysing filters 
can be easily designed. Wavelets derived by this method have potential application in the fields of optics, microwaves 
and electromagnetism. 
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1 Introduction 

A family of linear second-order differential equation with 
periodic coefficients was introduced by Emile Mathieu when 
studying vibrations in an elliptic drum 1 . Mathieu’s equa¬ 
tion is related to the wave equation for the elliptic cylin¬ 
der. This paper is concerned with the canonical form of the 
Mathieu Equation: 


antennas with arbitrary eccentricity, and (v) scattering by 
a coated strip. In general, the solutions of |l]) are not pe¬ 
riodic. However, for a given q , periodic solutions exist for 
infinitely many special values (eigenvalues) of a. For many 
physical solutions y must be periodic of period ir or 2n. 
It is also convenient to distinguish even and odd periodic 
solutions, which are termed Mathieu functions of first kind. 


i2 

d y 


^ + [a — 2q cos (2w) ] j/ = 0. 


(1) 


for a £ R, q £ C. The solution of Q is the elliptic 
cylindrical harmonic, known as Mathieu functions. Math¬ 
ieu wavelets as well as Mathieu functions can be appealing 
in a lot of Physics issues [2] including vibrations in an el¬ 
liptic drum, diffraction, amplitude distortion, the inverted 
pendulum, the radio frequency quadrupole, stability of a 
floating body, alternating gradient focusing, vibration in a 
medium with modulated density, and even when examin¬ 
ing molecular dynamics of charged particles in electromag¬ 
netic traps [3|[4]. They have also long been applied on a 
broad scope of waveguide problems involving elliptical ge¬ 
ometry [5 13 , including: (i) analysis for weak guiding for 


step index elliptical core optical fibres, (ii) power transport 
of elliptical waveguides, (iii) evaluating radiated waves of 
elliptical horn antennas, (iv) elliptical annular microstrip 
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2 Preliminaries and Background 


One of four simpler types can be considered: Periodic so¬ 
lution (n or 27r) symmetry (even or odd). For q ^ 0, the 
only periodic solution y corresponding to any characteristic 
value a = a u (q) or a = b„(q) has the following notation: 


Even periodic solution: 

c e u (w,q) = y A„ tTn cos (mw), for a = a v {q), 

m 

(2a) 

Odd periodic solution: 

se„(w, q) = ~y A^ m sin(mu>), for a = b v (q ), 

m 

(2b) 


where the sums are taken over even (respectively odd) val¬ 
ues of m if the period of y is ir (respectively 27 t). Given 
u, we denote henceforth by A m , for short. Elliptic 
cosine and elliptic sine functions are represented by ce and 
se, respectively. Interesting relationships are found when 
g-tcwo M): 

lim ce„(ro, q) = cos (uw), (3) 

g—5-0 

lim seRie, q) = sin(^w). (4) 

q —>0 

As it happens with trigonometric functions, Mathieu 
functions hold orthogonality properties: 
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The second non-periodic solution 14 


corresponding to 
c e u (w,q) is the Mathieu function Zce u {w,q). Similarly, 
The second non-periodic solution corresponding to s e v (w, q ) 
is the Mathieu function Z se u (w, q). 

One of the most powerful results of Mathieu’s functions 
is the Floquet’s Theorem [l5][l6|. It states that periodic 
solutions of 0 for any pair (a, q) can be expressed in the 
form: 


y{w) = F„(w) = e-> vw ■ P(w), or (8) 

y{w) = F v (-w) = e-i™ ■ P{-w), (9) 


0 1 2 3 4 5 6 


(a) 



where v is a constant depending on a and q and P(-) is ir- 
periodic in w. The constant v is called the characteristic ex¬ 
ponent. If v is an integer, then F v (w) and F v (—uu) are linear 
dependent solutions. Furthermore, y(w + kn) = eF kv y{w) 
or y(w + kn) = e~F k ^y(w), for the solution F v {w) or 
F v {—w), respectively. We assume that the pair (a, q) is such 
that | cosh(jV7r)| < 1 so that the solution y(w) is bounded 
on the real axis [T7| . 

The general solution of Mathieu’s equation (q £ R, v 
non-integer) has the form 

y(w) = c\e' vw ■ P(w ) + c 2 e~ :jL ' w ■ P(—w), (10) 


where Ci and c 2 are arbitrary constants. 

All bounded solutions -those of fractional as well as in¬ 
tegral order- are described by an infinite series of harmonic 
oscillations whose amplitudes decrease with increasing fre¬ 
quency. In the wavelet framework we are basically con¬ 
cerned with even solutions of period 2n. In such cases there 


exist recurrence relations among the coefficients 14 


(a-l-q)A 1 -qA 3 = 0, (11) 

(n to )A rn q(A m —2 T A rn j r2 ) = 0, (12) 


(b) 

Figure 1: Some plots of 27r-periodic 1 st kind even Mathieu 
functions. Elliptic cosines shape for the following set of 
parameters: (a) v = 1 and q = 5, (b) v = 5 and q = 5. 


as follows: 


x = acosh(u) cos(n), 

(13) 

y = asinh(u) sin(n), 

(14) 

z = z, 

(15) 


where u £ [0, oo), v £ [0,27r) and z £ (—oo,+oo). 

A further important property of Mathieu’s functions is 
the orthogonality. If a(v + 27t, q) and a(u + 2s, q) are simple 
roots of cos (f v) — y(f), then 


F v + 2 P (w) ■ F v+2s (-w)dw = 0, p^s. (16) 

In other words, (F„ +2p (w), F„ +2s {—w)) = 0, p / s, where 
(•, •) denotes the inner product. 
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for to > 3, to odd. 

Different methods for computing Mathieu functions are 
available in the literature [I8|[l9|. Figure [I] shows two illus¬ 
trative waveforms of elliptic cosines, whose shape strongly 
depends on the parameters v and q. 

3 Elliptic Cylindrical Coordinates 

The v coordinates are asymptotic angles of confocal hyper¬ 
bolic cylinders, which are symmetric with respect to the x 
axis. The u coordinates are confocal elliptic cylinders cen¬ 
tred in the origin. The transformation equations are given 


4 Multiresolution Analysis Filters and Math¬ 
ieu’s Equation 


Wavelet analysis has quickly developed over the 
years 


past 

20 and there has been a explosion of wavelet applica¬ 


tions including seismic geology, quantum physics, medicine, 
image processing (e.g., video data compression, reconstruc¬ 
tion of high resolution images), fractals, computer graph¬ 
ics, linear system modelling, computer and human vision, 
denoising, filter banks, radar, wideband spreading, turbu¬ 
lence, statistics, volumetric visualisation, metallurgy, solu¬ 
tion of partial differential equations, and so on. Essentially, 
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the wavelet transform is a signal decomposition onto a set 
of basis functions, which is derived from a single prototype 
wavelet by scaling (dilatations and contractions) as well as 
translations (shifts). 

In the sequel, wavelets are denoted by and scaling 
functions by with corresponding spectra U/(u') and 

&{w), respectively. The equation h n <j)(2t — n)., 

which is known as the dilation or refinement equation , is 
the chief relation determining a Multiresolution Analysis 


(MRA) 21 22 


5 Two Scale Relation of Scaling Function and 
Wavelet 


Defining the spectrum of the smoothing filter {hk} by 

(17) 

v z ke z 


the central equations (in the frequency domain) of a mul¬ 
tiresolution analysis are 


22 


*(„)=* (!).*(!), as) 

(19 > 


where 


G(w) = -7=^2 9ke 

V ^ L.c'Z 


—jwk 


( 20 ) 


is the transfer function of the detail filter. The orthogonal¬ 
ity condition corresponds to 20 22 ]: 


H( 0) = 1 and H{tt) = 0, 

( 21 ) 

\H(w) 2 + \H(w + tt) 2 = 1, 

( 22 ) 

H(w) = —e~i w G(w + tt). 

(23) 


6 Filters of the Mathieu MRA 


The subtle liaison between Mathieu’s theory and wavelets 
was found by observing that the classical relationship 20 


23 


*(«0 = e~^H* (!-.)* (|) (24) 


presents a remarkable similarity to a Floquet’s solution of 
a Mathieu’s equation, since H(w) is a periodic function. 
As a first attempt, the relationship between the wavelet 
spectrum and the scaling function was put in the form: 


'k(w) 

*(f) 



(25) 


Here, on the second member, neither v is an integer nor 
H(-) has a period tt. By an appropriate scaling of this 
equation, we can rewrite it according to: 


T(4-u;) 
$(2 w) 


= e~ j2w .H{2w 


TT). 


(26) 


It follows from © that H is periodic. We recognise 
that the function Y(w) = (4w)/H>(2 tx;) has a nice in¬ 

terpretation in the wavelet framework. First, we recall 
that 'k(w) = G( f)$(f) so that T(2u>) = G(w)$(u;). 
Therefore the function related to Mathieu’s equation is 
exactly Y(w) = G(2w). Introducing a new variable z, 
which is defined according to 2 z = 2w — tt , it follows that 
—y(z+§) = e~^ 2z H* (2z). The characteristic exponent 
can be adjusted to a particular value v as follows: 

- e-A "- 2 ) Z Y (z + |) = e~ jvz .H(2z). (27) 

Defining now P{—z) = H*(2z) = J2ke% c 2 k e ~^ 2kz -, where 
C 2 k — ^75 we figure out that the right-side of the above 
equation represents a Floquet’s solution of some Mathieu 
differential equation. The function P(-) is 7 r-periodic verify¬ 
ing the initial condition P(0) = ^ hk = 1) as expected. 

The filter coefficients are all assumed to be real. Therefore, 
there exist a set of parameters (aG,qc) such that the aux¬ 
iliary function 


y^) = -e- iiv - a)z Y v (z+^ 


(28) 


is a solution of the following Mathieu equation: 

d 2 y v 


dz 2 


+ (og - 2 q G cos( 2 z))y v = 0 , 


(29) 


subject to y„( 0) = —Y{tt/2) = —G{tt) = — 1 and cos(iz 7 r) — 
y v { 7 r) = 0 , that is, y v (-it) = (~ 1 ) 1/ . 

In order to investigate a suitable solution of (29), bound¬ 
ary conditions are established for predetermined a, q. It 
turns out that when v is zero or an integer, a belongs to the 
set of characteristic values a„(q). The even (27r-periodic) 
solution of such an equation is given by: 


2 b{z) = - 


ce „{z,q) 
ce„( 0 ,g)' 


(30) 


The Y v {w) function associated with y„(z) and related to 
the detail filter of a ‘Mathieu MRA’ is thus: 

Y v (w) = G„( 2w) = _ e A- 2 )(--f).^_M. (31) 

ce„ fu, q) 

Finally, the transfer function of the detail filter of a Math¬ 
ieu wavelet is 


G u (w) = —e 


__#(— 2 )(^) 


ce„( 0 , q) 


(32) 


The characteristic exponent v should be chosen so as to 
guarantee suitable initial conditions, i.e., G„(0) = 0 and 
G„( 7 t) = 1, which are compatible with wavelet filter require¬ 
ments [20p3] . Therefore, v must be odd. It is interesting to 
remark that the magnitude of the above transfer function 
corresponds exactly to the modulus of an elliptic sine 17 1 : 

|G„(w)| = se„ /ce„(0,q) . (33) 
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The solution for the smoothing filter H{-) can be found 
out via QMF conditions 21 , yielding: 


cefoO, q) 

In this case, we find H u (ir) = 0 and 

ce„(f ,q) 


\ H A W )\ = 


ce y (0 ,q) 


(34) 


(35) 


Given q , the even first-kind Mathieu function with char¬ 
acteristic exponent v is given by: 


ce v (w,q) = E A 2 i+\ cos(2 l + 1 )w, 


(36) 


;=o 


in which ce y (0, q) = ^2Z+i- The G and H filter coeffi¬ 

cients of a Mathieu MRA can be expressed in terms of the 
values {A 2 i+i}/ e z of the Mathieu function as: 


fef kl | 2 ;_ H /2 

V2 ce v (0,q) ’ 

9i = / 1 y ^|2i+^-2|/2 
y/2 ce„(0,g) 


(37) 

(38) 


It is straightforward to show that h u _i = h\ + and \/l > 0. 
The normalising conditions are ^ J2k=-oo = ~ 1 and 

Y^t=-oo(~ l) fc ^fc = 0- Illustrative examples of filter transfer 
functions for a Mathieu MRA are shown in Figure [2j for 
v = 1,5, and a particular value of q (numerical solution 
obtained by 5-order Runge-Kutta method). The value of 
a is adjusted to an eigenvalue in each case, leading to a 
periodic solution. 

Such solutions present a number of v zeroes in the inter¬ 
val |io| < 7 t. We observe lowpass behaviour (for the filter 
H) and lrighpass behaviour (for the filter G), as expected. 
Mathieu wavelets can be derived from the lowpass recon¬ 
struction filter by an iterative procedure (the same approach 
as the one usually used for plotting Daubechies wavelets). 
Infinite Impulse Response filters should be applied since 
Mathieu wavelet has no compact support. However a Fi¬ 
nite Impulse Response approximation can be generated by 
discarding negligible filter coefficients, say less than 10 -10 . 

In Figure [3j emerging pattern that progressively looks 
like the wavelet shape is shown. Waveforms were derived 
using the MATLAB wavelet toolbox. A fractal behaviour, 
which is common for some wavelets, can be noticed in these 
figures. As with many wavelets there is no nice analytical 
formula for describing Mathieu wavelets. Depending on the 
parameters a and q some waveforms (e.g., Figure. [3^,) can 
present a somewhat unusual shape. 


7 Conclusions 

A new and wide family of elliptic cylindrical wavelets was 
introduced. It was shown that the transfer functions of the 
corresponding multiresolution filters are related to Math¬ 
ieu equation solutions. The magnitude of the detail and 



Figure 2: Magnitude of the transfer function for Mathieu 
multiresolution analysis: smoothing filter \Hv(w) \ for a few 
Mathieu parameters, (a) v = 1, q = 5, a = 1.858; (b) v = 5, 
q = 5 a = 25.550; (c) v = 15, q = 5, a = 225.056. 


smoothing filters corresponds to first-kind Mathieu func¬ 
tions with odd characteristic exponent. The number of 
zeroes of the highpass |G(u>)| and lowpass \H(w)\ filters 
within the interval |u>| < n can be appropriately designed 
by choosing the characteristic exponent. This seems to 
be the first connection found between Mathieu equations 
and wavelet theory. It opens new perspectives on link¬ 
ing wavelets and solutions of other differential equations 
(e.g. Associated Legendre functions, Coulomb wave func¬ 
tion, Parabolic cylindrical functions etc.) Further general¬ 
isations such as the cases q < 0 or complex characteristic 
exponent can provide new interesting ’waves’. This new 
family of wavelets can particularly be an interesting tool 
for analysing optical fibres due to its symmetry |7 10 24 


Mathieu wavelets could as well be beneficial when examin¬ 
ing molecular dynamics of charged particles in electromag¬ 
netic traps such as Paul trap or the mirror trap for neutral 
particles |3l|4l- 
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Figure 3: FIR-Based Approximation of Mathieu Wavelets 
as the number of iteration increases: 2, 4, and 6 iterations, 
from the top to the bottom row. Filter coefficients holding 
\h\ < lO -10 were thrown away (20 retained coefficients per 
filter in both cases), (a) Mathieu Wavelet with v = 1 and 
q = 5 and (b) Mathieu Wavelet with v = 5 and q = 5. 
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